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LIST  CF  SYMZ30L3  IN  ORDER  OF  APFEARAVCE  IN  THE  TK>.T 


area  cf  one  elecrrode 

♦  Tickr.e?s  of  rho  elemort. 

polarlii.a^  ion  of  the  elt-ment 

A  P  =  inllia'i  charge  on  the  electrodes 
o 

vol tage 

electric  field 

Index  for  iruia’  ctaie  tne  shcck  front) 

index  for  final  Ftote  (ben.  tne  shock  froi.t) 
time 

=  normal  puise  length 
velocity  cf  the  shock  front 
current  output 

resistive  lead  (chmlc  resistor) 

inductive  load  (Inductance  of  the  resistor  and  the  wires) 
u  t  -  position  of  the  shock  frent  at  the  -ime  .* 
ciiarge  density  or.  tne  eleotrtdes 

C  A  “  total  charge  on  the  electrodes  ( as  a  furcricr.  of  time) 

Blec-ric  displacement  in  the  practical  systea. 

dielectric  constant  of  free  space 

dielectr.c  constant  behind  the  shcck  fren’ 

dielectric  constant  before  the  shock  front  (mean  value) 

»'  /I  ^  -  1  =  shape  factor 
if 

dielectric  constant  Defore  the  shock  front,  for  lew  field  E, 

coefficient  for  the  quadravic  fit  to  the  hysterejis  curve 

dielectric  constant  before  the  shock  front  (maximum  value  for 
coercive  foice  field,  E  -  -  E^) 

coercive  force  in  the  uncompressed  material 

t /t  -  normalized  time 
■  c 

I  t  /Q  =•  normalized  current 
o'  o 

(t  V/(c  AR) 

2' 

M(x  e  X  S  ''2) 
n/M/(2S; 
fSx> 
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T  =  (R  t  )/(2L) 
o 

N  =  (^^'t  2)/(e  lA) 

o  1 

Y  =  (1  -  x)/x 

(}l(x)=r  auxiliary  tunc ti mi 

C  =  capacitance 


w  ■=  en.'-rgy  output, 

r  =  total  lesistance  of  the  compressed  material 
p  =  resistivity  of  the  compressed  material 

(j  =■  charge  density,  diminished  by  conduction 

n  =  •  »-SV(.  ) 

2r  =  diameter  of  the  clement 
^  =•  position  along  this  diameter 

<p  =  angle  of  impact 

t^  =  time  delay  due  to  the  angle  of  impact 

z  =  t  /t  -  nc'rmalized  time  delay 

a  o 

B  =  abbreviation,  for  equation  (72) 

T  =  abbreviation  for  equation  (74) 

X  =  position  of  the  peak  of  the  output  curve  y(x) 
P 
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ABSTRACT 


When  a  polarized  ferroelectric  eiemert  as  traversedby  a  shock 
front  that  destroys  the  polarization,  a  current  output  appears.  The 
shape  of  *  he  curve  of  this  current  output  versus  time  is  computed  in 
a  iiormalized  and  general  form  for  *he  case  where  the  shock  front 
moves  an  a  direction  nooTnal  (or  arproxim i*:ely  normal}'  to  the  elec¬ 
trodes  and  parallel  to  the  p''larization  ''he  influence  of  the  external 
electric  lead,  of  the  dielectric  constants  and  the  hysteresis,  of  a 
conductiv’ty  in  the  compressed  material,  and  cf  an  oblique  impact  is 
discussed.  The  condition  for  maximum  en*'rgy  output  also  is  discussed. 
Mathematical  expressions  for  the  current  output  are  given  for  the 
various  conditions,  examples  are  calculated  and  sho»;  in 
Methods  are  described  for  evaluating  the  dielectric  ronctant  ir  the 
compressed  material  from  measured  output  curves, 

1 .  IN-RODUCTION 

Tne  calculations  reported  evolved  fr'^m  the  need  to  evaluate  cer¬ 
tain  experimental  results  The  topic,  experimental  setup,  and  first 
results  are  described  by  P.S.  Brody  1 n  ref c lence  1,  in  which  he  states 
the  basic  theoretical  expression  for  the  current  output  under  shock 
compression,  when  a  shock  front  travels  through  an  axially  polarized 
ferroelectric  disc  in  the  direction  of  i‘"i  polarization  [n  the 
present  report,  the  theoretical  expression  (transfonned  into  a  more 
convenient  and  general  ''normalized  form")  is  examined  in  terms  of  its 
various  pai'ameters.  New  parameters  are  added  --  the  curvature  of  the 
hysteresis  curve,  the  inductive  load,  the  conductivity  behind  the 
s.ncck  front,  and  the  angle  cf  impact.  Ir  certain  cases  the  formulas 
can  be  solved  for  Cf,  the  dielectric  cons  .ant  of  the  compressed  mater¬ 
ial.  The  latter,  therefore,  can  be  calculated  from  measured  output 
curves , 


The  assumptions  are: 

1)  The  shock  front  1?  very  thin,  tf'v  is:  pressure,  temperature, 
particle  velocity,  and  dielectric  constar,-.  reacn  their  new  values 
vithin  a  time  short  compared  with  the  tota'  pulse  length. 

2)  No  relaxation  time  wit.h  respect  to  the  change  in  polarization 
appears.  Tne  polarization  goes  to  zero  immediately  after  the  material 
is  in  its  high  pressure  state. 

3)  The  material  is  isotropic. 

4)  No  reflected  rarefaction  or  compr-'ssion  waves  from  the  edges 
or  from  the  back  surface  of  the  element  a,.pear  in  time  to  disturb  the 
process . 

5.  The  uncompressed  material  (before  the  ^hock  front)  xullows  a 
hysteresis  curve,  which  can  be  measured  fc.  elements  of  ihe  same  type, 
and  for  which  therefore  a  iiumerlcal  fit  can  be  made  for  the  purpose  of 
evaluating  measured  output  curves. 
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6)  The  "dielectric  constant"  of  the  compressed  material  (be¬ 
hind  the  shock  front)  is  supposed  to  be  a  constant,  with  respect 
to  the  field,  Ej .  If  this  assumption  is  not  true  for  a  large  field 
(which  we  do  not  know),  tlien  the  early  part  of  the  output  pulse  will 
be  different  from  that  shown  in  this  report.  The  tail  of  the  pulse, 
however,  def*cnds  rm  i  r  lor  a  small  lield,  when  it  is  more  likely  to 
be  constant. 

2 .  CURRENT  OUTPDT  AS  A  FUNCTION  OF  TIME 

In  order  to  .«et  up  ^in  equation  ‘or  the  current  output,  the  fol¬ 
lowing  physical  process  i.-<  do.*- c  ri  ted :  Initially  tiic*  wnole  element 
is  polarized,  the  value  of  the  polarization  is  P  .  A  free  charge 
±Q  (-•:  A'P  )  has  been  allowed  to  .iccumul  at  e  or.  tRe  dec  t '"odes ,  kept 
in  place  b/  the  polarization  ar.d  neutralizing  it?  c'loct  The  net 
field  is  zero  (fig.  la).  After  an  impact  ha.s  occurred,  a  shock  front 
travels  tnrough  the  element,  leaving  behind  a  cubit:  lattice  with  no 
spontaneous  polarization.  During  this  time  the  charge  from  the  elec¬ 
trodes  flews  off.  The  internal  polarization  in  the  undasturbed  region 
is  no  longer  completely  neutralized;  an  increasing  field  E  appears  here 
in  the  negative  direction.  In  the  compressed  region  with  no  spontan- 
eou.s  polar  ization ,  a  field  Ej  in  the  positive  direction  appears  as  a 
result  of  the  external  charge  on  the  electrodes.  Ej  decreasc-r  with 
time  as  the  charge  does  (fig.  lb).  Finally  the  shock  front  has  passed 
the  element.  Tnere  remains  a  decreasing  field  Ep  as  long  as  a  charge 
remains  on  the  electrodes  (fig  Ic) . 

Because  the  thickness  of  the  element  is  small  compared  with  its 
diameter,  and  because  the  dielectric  constants  in  the  material  are 
much  larger  than  in  free  space,  edge  effects  can  be  neglected.  The 
sum  oi  voltages  along  a  closed  loop  must  be  zerc.  So  we  get  (see 
fig.  lb): 

•  (f -C'  “  I  •  R  “  L  ^  -  0.  (1) 

E  is  negative.  The  current  I,  in  the  direction  .shown,  is  the  nega¬ 
tive  derivative  of  the  charge  Q  with  respect  tc  time  (since  the  charge 
decreases  with  time).  In  the  practical  system  the  charge  density  O 
equals  tne  dielectric  displacement , O  -  D.  D  is  conneoted  to  E  in  the 
foil' wing  ways;  For  the  compressed  material  we  assume  a  linear  rela¬ 
tionship,  employing  g p ,  the  dielectric  constant  behind  the  shock  front, 


D  -  gp  Ef  , 


(2) 


from  which  follows 
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Tiif  un('<mpri‘ssi'(i  material  l)efo,’o  the  shock  front  follows  ^  h^’^teresls 
iiirvc-.  We  are  interested  111  Its  second  quadrant  only. 


^  m 


For  a  very  small  field  E  we  find  the  relation* 

D  =  P  +  a  E, 
o  ' 

a  Is  the  apparent  dielectric  constant  when  measured  with  a  low  field 
bridge.  With  increasing  field  £  in  the  negative  direction,  the 
AD 

dielectric  constant  —  is  no  longer  constant  but  increases  also, 

and  reaches  its  maximum  value  e  for  the  field  -E„  (coercive  force). 

m  c 

In  a  simplified  manner  the  hysteresis  curve  can  be  approximated  by 
a  (somewhat  arbitrary)  mean  value  ,  giving 


*  The  difference  between  polarization  P  and  dielectric  displacements 
D  is  neglected  since  is  very  large. 
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€  E 
1  1 


D  =  Fo 


from  which  follows 

P  -D 

„  o 


1  e . 

1 


€  . 

1 


(3> 


A  better  approach  is  accomplished  with  a  second  order  fit, 

D^P  +aE.  -6  ,  (4> 

o  1  i  ’ 

from  which  follows 


The  current  is 


and  the  distance  traveled  by  the  shock  front  is  ^  u  t.  Combining 
(1)  with  (2)  and  (3)  results  In  ^ 


_Q 

^f 


u  t 
s 


u^t)  4-  ARQ  4  ALQ  =  0 


(5) 


This  is  a  differential  equation  of  charge  Q  as  a  function  of  time  t. 
From  its  solution  one  gets  the  current  I  =  -  Q  as  a  function  of  time. 

P  >  u^ ,  R,  A  are  constant  experimental  values  or  parameters. 

T^e  Initial  values  are  Q  =  Qq  for  t  =  0  and  Q  =  0  for  t  =  0. 


Since  this  report  deals  with  the  general  theory  and  the  influ¬ 
ence  of  certain  parameters  rather  than  with  specific  values  from  exper- 
imeiiLs,  it  is  convenient  to  eliminate  some  geometrical  and  experimental 
parameters  by  introducing  a  normalized  form.  Nothing  is  lost  with  it, 
since  any  specific  values  always  can  be  introduced  again  into  the  re¬ 
sulting  equations.  Normalization  in  this  sense  is  achieved  by  normal¬ 
ized  units  of  time,  charge,  and  current.  Tlie  unit  of  time  is 


t 

o 


(  = 


f 


the  time  required  for  a  shock  to  traverse  the  element.  Tne  unit  of 
charge  is  Q  ,  the  initial  charge  bound  to  the  electrodes  by  the  polar¬ 
ization  of  ?he  element.  The  unit  of  current  is  tlie  quotient  ~  .  So 

we  get  the  normalized  time 
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*lM(  1  1  rum  /iiii  ii>  mw  (ixctpl  iitidcr  an  oljl  li|Ui'  impact,  where 

"  <  X  '  1  *  /  )  ,  .111(1  111"  IK  inii.i !  1 /.<■('.  call  ii  lit  ..ui.'jul 

t  I. 

II  I)  clQ  1  cJU 

Q  fit  Q  dx 

Miisi  t;ra:>h.s  in  this  repoi  t  .-.how  l ‘ii-  fiiiitiKMi  y(x). 

Hfiuatlon  (5)  multiplied  by  gives 

Vj 

f  AR  'i.AL 

-Q  (1-x)  I  Q  (1+S  x)  +  Q  — n —  +  Q  — y-  -  0 

/O  L 


with  S  i-.  -i  -  1  ~  "shape  laeter". 


In  order  to  use  the  normalized  time  x  rather  than  t,  we  have 
to  t  ransform 

Q  ^  ,  y  „  _L  . 

^  t  dx  ^  2^2 

o  Iq  dx 

Riving 

£  AR  2  e  AL 

-«^(l-x)  »  ,  <1*8  ,)  ♦  jS  ^  ^  ^  0  .  ,6) 

®  dx  £ 

This  Is  the  general  differential  equation  for  Q(x);  the  initial 

do 

values  arc  Q  =  and  ^  =  0  lor  x  r;  0.  Once  the  solution  Q(x)  has 

been  found,  the  normalized  current  output  y(x)  Is  computed  by 

1  dQ 

For  the  case  R  =  n,  L  0,  {short  clriuit)  we  got 


Q  =  Q. 


o  1+Sx 


(l  +  f5x)‘ 


For  the  case  R  ^  0,  L  =  C  we  get  from  (8) 


^  +  0  (1  +  S  x)  M  - 
dx 


,1  (1  -  x)  M  ^  0 
o 
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A  raiher  simple  solution  exis's  if  S  -  0  (or  ^  ;;  ^f)>  namely 


«  -  Q„(> 


-Mx 


1  c 
^  M  ~  M 


(12- 


f rom  which  foil ows 


,  -Mx 
y  =  1  -  e  , 


(13) 


and  the  initial  slope  is 


I 


,  dx 


=  M  = 


^X=:0 


e.  AR 
1 


(14) 


If  S  ^  0,  the  general  solution  to  (11)  follows  the  scheme: 


with 

Here  we  have 


Q'  +  Q  f(x)  -  g(x)  =  0 

,  X 

Q  =  e  '  (Q  +1  g(x)  e*^  dx) 

Jo 

F  = I  f (x)  dx  . 

Jo 


S  2 

F-M(x+*x  ),  and  thereby 


0  =  <J„  e 


-F 


[*  ♦  “/' 

•^0 


(1  -  x)  e  dx 


]• 


(15) 


y  =  M  l^d+S  x)  e  ^  +  M(l-rSx)  e  (i-x)  dx  -1  +  xj  .(16) 

The  integral  cannot  be  solved  analytically;  however,  by  partial 
integration  another  form  can  be  derived: 


y  =  M 


i±§ 

s 


(1-fSx)  e  ^  -f  M(l4-Sx)  e 


Jo 


<16a> 


A  third  form  of  the  solution  can  be  derived  from  (16)  by  these 
definitions: 


0)  =  v-5~  and  to  =2  w  (1+Sx) 

O  aS  o 


from  which  follow 
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/  MS 

dco  -  Si*)  dx  -  :  --  d>. 

C>  V 


S  2^  2  2 

F  .  M(x  +  -  X  )  -  (1)  -  00 


1  +  Sx  = 


1  +s 

1 

(0 

"s 

S 

CO 

f> 

CO 

CO 

o 

y(x  )  =  y  (  u)(x))  y(co) 


y(oo)  M 


2  2 

J0_  ^-(00  -0^ 


00, 


1  JO 

s  <0^, 


2  2  /  \22,  -i 

M  -•‘-o  >  (its  _  i  .(0  ]  (0-<o„  _l  _ 

00  I  ^  S  S  CO  /  Soo 

O  j  o  o  _J 

00^ 


which  after  some  manipulations ,  leads  tc 


y(<o>  =  M  — g— 


,2  2.  2  r  “  2 

JO  ^-((O-Ot, 

00 

00 


00  .  , 

e  dco  '  1 


(leb) 


Solution  (i6b)  is  given  in  reference  1  Its  advantage  is  that  the 
2,  fo  2 

“(0  I  00  I  j  -I  • 

term  e  I  ^  caii  be  tound  tabulated,  at  least  in  a  certain 

*^0 

range  of  uj. 

The  initial  slope  can  be  calculated  frcm  any  of  tnc-si  forms,  giving 


=  M  (1  +  S)  =  ■ 


x=0 


f  jPA 


(17) 


Equation  (17)  includes  equation  tl4)  for  S  0  I’rans  fr  rmea  bark  to 
the  observable  output  curve  l(t),  we  have  as  initial  ‘•lofe  for 
any  S: 


/<a\  ^  i! 

'•'■vV^o  ■  ’  '< 


(  lb  ) 
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'  .>11  age 


It  IS  inteiesting  to  note  that  the  initial  rise  of  curve 
(across  R)  versus  time  is  independent  of  R, 


Vdt4  .  ~  ■  t  €,A 
t=o  of 


This 


is  all  for  the  case  R  ^  0,  L  =  0. 

For  .the  case  R  =  0,  L  ^  0  we  get  from  (8); 

^  +  N  (1  4  Sx)^Q  -  N  (1  -  x)  Q  e  0 
,  ^  *  o 

dx 


with 


<^t  " 
N  ^  o 


e.LA 

I 


A  solution  can  be  found  if  S  = 


Q  =  (I  -  X  +  — 


0,  namely 
sin  x) 


(19) 


(20) 


(21) 


y  =  1  -  cos  >/  N  X 


Equation  (23)  is  true  for  any  S,  since  it  is  a  consequence  of  (20). 
For  the  general  case  R  4=  0^  L  4^  0  we  get  from  (8): 


with 


and 


Solutions  can  be  found  if  S  =  0,  namely 


d^Q 


— f  +  2T  g  +  N  (1  +  Sx)  Q  -  N  (1  -x)  =  0 

dx 


Rt 


T  = 


2L 


N  = 


e^LA 


(22) 

(23) 


(24) 
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^  I  K“  -J  / .  2  , 

<1 )  it)r  1  <  N  ,  1 1  p  sj  N  -  T  _  Q 

I,  A  ' 

1 


.2  ^2. 


^  ~  -^u  '^n'  "  '•  "  ‘'  "“no 


~Tx  2T  A 
e  —  cos  (j>x 


y  1  -  C- 


T\  T 


-T\ 


^  •■'in  ijx  -  c  os  (j)x  , 


(25a) 

(26a) 


t)  tor 

L*  A 

1 


Q  =  Q, 


2+r  -Tx-  ,  2. 

'*  ♦  T>  . 


(25b) 


y  1  -  e  C  Ibi  +  1 ) 


(26b) 


2  4  P  i  2 

c)  for  T  >N,  ~  >  ~>  *ithV’r  -  N  = 

Mj  t  a 
1 


9 


=  Q  [~J— 

o  I  N 


The  -(T-e)x  T-e  '(T+e)x 
"  2e(T-e)  26(1+9) 


(25c) 


y  ^  1  -  e 


-Tx 


T+9  ex  T-e  -Ox 
29  ■  29 


(26c) 


The-  Initial  siojn  is,  as  a  consequence  of  equation  (24), 


and  therefore  is  true  tor  any  S. 

So  tar  the  calculations  have  involved  the  simplified  formula  (3) 
for  the  hysteresis  in  the*  uncompressed  material.  There  has  also  been 
given  a  better  approximation  to  the  curved  hysteresis  with  formula  (4). 
The  reasons  for  using  t^’”  simplified  formula  (3)  at  all,  is  that  with 
formula  (4)  the  mathematical  expressions  become  very  clumsy,  the  cal¬ 
culations  cannot  be  pushed  so  far,  and,  on  the  other  hand,  many  basic 
features  can  be  demonstrated  adequately  by  the  simpler  formulas. 
However,  it  is  necessary  to  see  how  the  formulas  and  the  current-output 
curves  actually  look  when  the  curvature  of  the  hysteresis  is  taken 
into  account  better.  Because  of  the  mathematical  difficulties,  this 
will  be  done  for  the  case  R  =  0,  L  =  0  only. 
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liistcail  ()1  cfiuiiiion  (5)  w<'  tjy  combi  n  i  ll^;  (1)  wl  (h  (2)  and  (4)- 


u  t 


U 

2fl 


7 


I  Afi 

1  +  ^  (P  -  a) 
2  o 
CC 


(  £  -  u  t)  =  0. 
s 


(28) 


Coinparod  with  equation  (5),  the  last  two  terms  involving  R  and  L  have 
l)ecn  di'oppeu  111  ready. 


Th.' 

I  V 1  tig 


rc  1  a  t  ■.  on 


u  t 
s  o 


is  used  and 


the  time  Is  normalized, 


t/t^, 


._a  X  - 


1  -a) 

a 


(1  -  x) 


(29) 


For  X  =0,  we  have  =  P  and  therefore  from  (4);  E.  =  0,  and  from  (2) 
’  o  i 

P 


For  X  Vf"  0,  (29)  is  rearranged  to 


o  +  <T  - 


€jY  (Cl  + 


=  P 


2  2 
€,  Y 


with  the  solution 


€^Y  (CH-e,Y) 


23 


43P 


(QJfe^Y) 


(3Q) 


l*x 

including  the  abbreviation  Y=  - 

^  X 


Equation  (30)  gives  the  charge  per  unit  area  remaining  on  the 
electrodes  as  a  function  of  time  ft,  or  x  =  t/t^,  or  Y  =  ). 

For  the  very  first  part  of  the  pulse,  that  is  for  x  very  small 
and  Y~i-*oo,  a  simplified  formula  instead  of  (30)  can  be  deduced.  The 

X 

root  is  expanded  according  to  the  scheme 

2 


vTT 


1  z  z 
^  2  "  8 


giving 


C=  P 


-^-p2 
g  o 
CfY 


2  2,,  a 
CfY  (U^ 


J 


'll  U'li  111  turn  IS  (‘xpandod  again.  Dropping  all  terms  with  !iclo''s  o( 
2 

r  grca'er  than  x  ,  ve  obtain 

a  2  Cl  2 

o  r  -  p  --  X  -  p  X 

I  >  <  >  I  o  ^ 


'  'I '  I  >111(1 1 1  >  i‘  so  1 11 1  I  <  in  'tills  is: 


Q  i  P  ,  I'or  X  =  0 
o 

«2ft2 

a  -  P  -  P  —  X  ~  P  — X  for  X  very  small 
o  o  e  o  2  ’ 


rjYCai-fjY) 


O  -  0,  for  X  =;  1 


(OH-e^Y) 


~  1  »<■''*  /V  <  X 


Now  tlie  fields  can  be  calculated.  Using  equation  (2)  and  (31)  we 


E,  =  —  ’  for  X  =  0 


E  -  P  — ^  X  ,  for  X  very  small 

t  >.■  o  2  ’ 


Y(a+e  y) 

"f  - 


Ej.  =  0  ,  for  X  =  1 


-  1  ,  for  0  <  X  <  1 


(a+€jY)  j 


E.  Is  defined  by  equation  (4).  We  take  advantage  of  equation  (29)_, 
Which  Is  nothing  but  E^x  +  E^  (1  -  x)  =  0,  or  ~  l/Y  Ej.. 


f  K  -  n 


,  for  X  =  0 


E,  =  -  —  X,  for  X  very  small 
'  'f 


CM-;.r 

^1  =  -  -2— 


E,  =  - 


(ote^Y) 


-  1  ,  for  0  <  X  <  1  (33) 


wa  th  Y 


1  +  — --  -  1  ,  for  X  r.  1  , 


♦  Cf=  0  and  E.  0  for  x  =  1  is  typical  for  the  short  circuit  case 

(R  =  L  =  0) 
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'!  Ill'  (  ii;i  r;,!'  i  s 
ion  i.'i  I  I  /.(  (I  cii  ri'i'ii  I  I  ■; 


AO,  and  i  hi' i  c  I  orr  is  civoii  by  (  d  I  >  The 


V(x  ) 


J_  'JS 

Q  (lx 


1  (1  a 

P’  (Tit 

o 


1_ 

P  X 


_ 

2  iJY 


o 


» y  i  r-  Mil  fijua  r  i  mi  (Ml): 


y  ^  —  ,  for  X  =  0 
‘  ( 


—  f  P  X  ,  for  X  very  small 

f  °  ‘■f 


y  = 


26P  x^  L 
o 


with  yj 


'(—  +y)('/'  -  1)  -  Yd  -  Jr  ) 


,  for  0  <  V  <  1 


43P 


1  +  - ^ - -  and  Y=  — ^ 

.2  X 


(q:+€jY) 


cte. 


y  = 


23P 


1  +  - 1  ■  ^  >  for  X  =  1 

a  / 


=  —  ,  for  X  =  1  (as  is  explained  below), 


For  X  =  1,  equation  (31)  gives  0=  D  =  0,  therefore  we  must  con¬ 
clude  from  figure  2  that  E.  =  -  E  for  x  =  1.  By  (33)  this  is 

1  c  _ 


r,  „  Ct 

E  =  -  E  =  -  — 

i  c  23 


43P 


1  + 


-  1 


a 


Again  from  figure  2  we  find 


a 


e  "  P  ~  23P 
m  o  o 


43P 

1  +  -2-  -1 

a 


which  leads  to  yd)  =•  — 

€ 

m 


(34) 
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Again,  as  in  section  1 ,  we  start  with  the  si.aple  assumption  (3)  rather 
than  (4).  The  first  case  to  be  considered  is  that  of  a  short  circuit 
R  =  L  =  0: 


It-S 

( 1+Sx)^ 


(10) 


Here  the  only  parameter  is  S  =  -  1.  Figure  3  shows  a  set  of  curves 

for  different  3.  For  S  -  0,  the  current  output  is  constant.  Our  exper- 
'raents  gave  curves  of  thtj  type  S  >0  (i.e,  Cj  >  e^)  only. 

For  the  case  R  ^  0,  L  =  0,  we  have  the  solution  (13)  for  S  =  0 
and  (16)  for  S  ^  0.  Figures  4  and  5  show  cui^es  for  S  =  0  and  S  -  1, 
respectively;  the  paiaiucter  R  has  the  values  0,  10,^100,  1000,  10,000, 
100,000  ohms.  The  limiting  value  R  =  0  is  described  by  equation  (10). 
For  the  other  limiting  value  R-*’a)  (for  any  S)  wf  get 

■it 

y(x)  =  M  (1  +  S)  X  = 
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■  in'-  iiic-aiis  lhat  fir  an  open  circiili  the  curve  y(>>  for  i  rcreasi  t.rr  x 
t;oe£>  on  rising  wi  in  li.s  Initial  slope,  as  given  by  equation  ('17', 

Kqu.iiion  (35)  can  be  found  In  the  following  way  (using  M  =■  0  for 
ii  ‘ii>).  (f  S  =  0,  we  get  from  (13)  y«  1  -  (1  Mx>  =  ibr.  If  S  4  0. 
we  t.'il;.-  equaMnn  (iPn'  ,  put  !n  I'  M  (x  <•  ®  x^)  =  0  and  ge* 

V  =  M  (1  t  S)x  ^ 

It  should  be  mentioned  lhat  the  curves  <.>f  figure  5  could  have 
beef  plotted  'sing  one  of  the  equations  (16),  a  tedious  task.  DOFL's 
analog  computer  was  used  Instead,  which  employs  the  dl f feren tl .•1  .-qua- 
'lon  ( 11 ;  rather  than  Its  solution. 

The  next  case  to  be  discussed  Is  that  for  R  4=  0,  L  4  0-  ^ 

the  fixed  value  10  ohms  is  used  (wnich  corresponds  With  ihe  experiiuentsl  . 
Figure  6,  with  S  =  0,  shows  a  curve  for  L  =  10*'^  henry  and  Is  plotted 
according  to  equation  (26a).  Figure  7,  with  S  =  1,  shows  curves  for 
h  =  0  5  lO'^  and  1  -  10~  henry.  They  are  plotted  by  IXlFL's  analog 
computer,  using  the  differential  equation  (24).  With  all  these  curves 
we  note  oscillations  -vhlch  are  damped  after  one  or  two  periods  super¬ 
imposed  on  the  basic  curves  for  R  =  10  ohms  and  L  =  0.  Higher  values  of 
L  (not  s.hown  here)  result  in  oscillations  of  still  higher  amplitude  and 
less  dampliig,  until  eventually  the  whole  curve  shape  is  masked  by  these 
oscl  1 1  aiiofis .  As  far  as  the  experiments  are  concerned,  the  inductance 
1.S  a  disturbance  and  should  be  as  small  as  possible.  The  values  for  L 
shown  in  figures  6  and  7  are  about  those  that  enn  be  expected.  Th^can 
infliienoe  only  the  very  first  part  of  the  pulse. 

Introducing  the  quadratic  fit  (4)  for  the  hysteresis  of  the  un¬ 
compressed  material.  Instead  of  the  linear  fit  (3)  as  used  hitherto  in 
this  section,  we  discuss  the  short  circuit  case  R  =  L  =  0.  In  order  to 
point  out  the  difference  between  (4)  and  (3),  the  following  results  have 
to  be  compared  with  the  former  results  with  equation  (10)  (see  fig.  3). 

The  calculation  employing  the  quadratic  fit  (4)  Is  in  better  agreement 
with  experiments  (ref  2).  Using  eaution  (34),  three  euWes  have  been 
calculated  a::d  plotted  In  figures  8,  9,  and  10.  The  values  of  a.  Pg 
(as  given  In  the  beginning  of  this  section)  were  measured  statically  with 
elements  of  the  type  used  In  the  shock  experiments.  A  dynamic  measure¬ 
ment  of  a  and  3,  Imitating  the  actual  shock  conditions.  Is  planned.  So 
far,  however,  no  better  values  are  known.  As  dielectric  constants  behind 
the  shock  front  have  been  assumed;  Cj  =  1250  -  £„  (jflg.  8),  =  1540  ■ 

(fig.  9),  and  ej  =  2000  •  (fig.  10).  Using  a  mean  value  (according  to 
(3))  of  =  4400  •  this  corresponds  to  S  =  2.5,  1.9,  and  1.2  respec¬ 
tively  (see  fig.  3  for  comparison). 

With  equation  (10)  we  had  the  following  situation:  The  curve  starts 
with  y(0)  i  Cj/ef  and  decreases  steadily  with  increasing  x  until  y(l>  = 
£f/£^  .  With  equation  (34)  In  essence  the  same  thing  hanoens,  except  that 
Is  not  a  constant  (see  fig.  2). 
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R  =  10  ohms)  . 


No#  e.  starts  with  the  value  a  and  ends  with  the  value  e  .  Therefore 
the  cuWe  starts  with  the  y<0)  = 

between  the  curve  shows  a  maximum,  ilie  curve  also  he.s  a  marked  bend  in 
its  rise,  which  in  some  cases  could  be  observed  in  the  experimental  re¬ 
sults,  but  in  most  cases  was  masked  by  other  effects.  Die  meaning  of  the 
auxiliary  function  (k(x)  will  be  explained  later. 


4.  CURRENT-OUTPUT  BEYOND  t  s  t^ 


So  far  .'othlng  has  been  said  about  the  current  output  for  t  >  t,,, 

l.e,  for  the  time  after  the  shock  front  has  reached  the  back  surface 

of  the  element.  The  free  charge  on  the  electrodes  at  the  time  x  »  1 

(t  =  t  )  is  0(1)  and  can  be  calculated  from  the  appropriate  equations, 
o 


If  there  is  no  disturbance  by  further  shock  or  rarefaction  waves 
or  reflections  from  the  side  walls,  and  if  th3  connecting  wires  dc  not 
break,  the  element  then  is  nothing  but  a  capacitor.  It  has  a  unifoni 
(compressed)  material  with  the  dielectric  constant  €.f.  Neglecting  the 
slight  geometrical  compression  (shrinkage  of€.),  the  capacitance  is 


Remembering 


and  1  -f  S 


H  (1  +  g> 


we  find  the  identity 


which  is  the  normalized  decay  time  of  the  element  as  a  capacitor.  For 
the  time  x  >  1  this  capacitor  discharges  according  to  the  law 


wx)  .  q<i)  •  «  S)  <»  - 

With  y  s  -  —  the  noraialized  discharging  current  is 

Qo  ‘•x 

y  ^  forx>l 


(36) 


(37) 


The  total  charge  outp'it  is 


CO  ,oo  ( 


dQ  «  I 


(38) 


which  is  equivalent  to 


(39) 
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i'(ir  sditic"  appl  icat  i(jns  the  quest  ion  arises  as  to  which  value  of  R 
ti'suMs  in  a  maximum  of  enertjy  output,  W. 


0 


this  prohiein  can  be  calculated  analytically  fc'r  the  case  S  -  0.  The 
(iofinitioii  of  the  effective  pulse  length  T  is  somewhat  arbitrary,  and 
tie  pends  on  the  experimental  conditions.  We  put  T  =  t  ,  because  the  ex  • 
(lerimoMs  infiicule  that  the  discharge  virtually  comes  ?o  an  end  at 
due  If  ra  n  f  ac  t  1  I'll  waves  and  subsequent  breaking  c'l  i  he  mres.  i. 
must  be  noted  that  in  general  not  all  the  free  charge  from  the  element 
can  flow  off  during  this  time. 

We  limit  the  calculation  to  the  case  R  0 ,  L  -c  0,  S  =  0,  with 
•he  current  output 


% 


y  =  1  -  (' 


•MX 


wl  th  M  - 


t  .  AR 
1 


from  which  foilows 


Q  Q  -  -^-1 

I(t)  =  *-py»--p(l-e  °  ) 


W  =  J  I  h  dt 
0 


w  ^ 


c  .  A 
1 


”  '  2m2  *  " 


2  -M  1 

2 


2M 


(13) 


(53) 

(52) 


(54) 


W  becomes  a  maximum  for  a  certain  value  of  R  cr  M,  which  is  determined 
dW 

by  putting  ^  =  0-  This  leads  to  the  equation 


*--M  '2M 

3  -  M  -  (4  +  2M>  e  +  (liM)  e 


0 


with  the  solution  M  =  1.8926.  With  this  value  the  maximum  energy  output 
is,  according  to  (54) 

Q 

^max  a  ~~~7  ■  0.19  (55) 

e.A 
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Ihis  is  only  38  percent  of  the  pc>tential  energy  of  the  element,  if  it 
Is  considered  a  capacitor  of  capacitance 


C  = 


which  has  the  potential  energy 


W 


Q  ^  Q 

Jo _ JOJJ 

2C  “  e^A 


•0.5  . 


The  total  charge  output  up  to  the  time  t^  is  according  to  equation  (43) 


«o  -  -  5  •  i  =  <io  •  <“> 

Using  the  specific  values  from  oiir  experiments  as  introduced  in  section  2, 
1264 

we  get  M  =  — —  which  in  '•omhi nation  with  M  =  1.802C  leads  to 
R  =  668  ohms 


as  the  value  of  the  external  resistor  for  which  the  energy  output  is  a  max¬ 
imum.  The  latter  is  according  to  (55) 

\ax  =  Joules. 

6 .  (XINDUCTIVITY  BEHIND  TWS  SHOCK  FRONT 


So  far  the  assumption  has  been  made  that  a  ferroelectric  material  is 
nonconducting  behind  as  well  as  in  front  of  the  shock  front,  at  least  for 
a  pressure  below  125  kl lobar  and  temperatures  below  ISO^C.  However,  there 
is  at  least  one  reference  in  the  literature  (ref  3)  stating  a  specific 
resistivity  p  s  1000  ohm  *  cm  for  a  pressure  of  300  kilobar  in  barium 
titanate.  As  the  following  calculations  will  show,  a  resistivity  larger 
than  p  *  10®  ohm  •  cm  should  not  be  detectable  in  the  output  curves  (in  a 
short-circuit  situation),  and  nay  therefore  be  regarded  as  infinite.  The 
influence  of  a  lower  resistivity  is  the  topic  of  the  following  calculations. 
The  model  used  in  section  2  (fig.  1)  is  introduced  again,  with  L  a  O  and 
R  very  small  (short-circuit  assumption).  The  linear  fit  to  the  hysteresis 
curve,  equation  (3),  is  used.  The  conductivity  behind  the  shock  front  is 
treated  as  a  shunt  resistor  with  the  resistance 


r  ap  e 
A 


9 


P  is  the  specific  resistivity  of  the  material  under  compression. 
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Figure  11.  Element  as  the  shock  front  travels  through  It,  and 
equivalent  circuit. 

In  the  former  treatment  (section  2)  it  was  supposed  that  no  free 
charge  can  accumulate  at  the  interface  between  the  two  regions  of  the 
element  (in  other  words:  at  the  shock  front),  which  Is  equivalent  to 
the  statement  that  the  dielectric  displacement  D  (equal  to  the  charge 
density  a)  is  continuous  through  this  Interface. 

Now  we  have  a  shunt  resistor  which  draws  a  certain  current  i, 
therefore  D  is  different  in  each  region  of  the  element.  is  a  in 
the  region  before  the  shock  front.  Behind  it,  Df  is  wl.ich  is 
smaller  than  a  . 


The  basic  equation  is  Uj  +  U2  +  U3  =  0,  or  with  short-circuit 
assumption  +  U2  =  0.  We  look  at  the  regions  of  the  element  as 
capacitors,  with  capacitance 


It  follows 


p  S  eA 

^  "  U  “  U  "  T 


U  =°d 


o-P 


U.  = 

JL 


a-P.  £ 


d-o  =  — (t  -  t) 

^  €4  t  O 

1  1  o 


and  also 


U„  =  -^5=-^^t 


2  “ 

I  ^  O 

4 

Ug  =  IT  =  i^ 


3'. 


! 


Till!  d  1  spl  ac’omeiu  cur  rent  Is 


acrijss  the  front  region  1  -  A  — "  and 

da 

ficnisc:  the  back  region  I-i  =  A  — , 
from  wliicli  it  follows  that 

dP  dO, 

Tlie  two  expressions  for  combined  result  In  a  differential 


equation  ford  ; 


da  1  _  d<7 

dt  ^  pCf  ^  ~  dt  ’ 


whose  solution  through  the  initial  value  a  =*  a  ^or  t  s  0  is 


0=0 - e 

P€f 


— 

oef  [ 


a  e  dt 


and  can  be  used  to  eliminate  o. 


Introduction  of  this  a  into  any  of  the  two  equations  for  Ug  leads  to 


u  _ il- 

2  €ft-  2  ^ 

‘^fPto 


t  t 

P€f  f  ^  P«f 
o  e 

0 


Now  the  basic  equation  is 


U,  +  U„  =  - -  (t  -t)  +  -°-  -r^  t - ^  ^ 

^  2  t^  o  ej  t^  ^  t^ 


t  t 
P^f  f 


t  t 

Wi 

a e  dtsO 


whlcli  must  be  solved  fura(t).  It  may  be  noted,  by  comparison  with  equation 
(1),  that  the  third  term  only  expresses  the  disturbance  due  to  a  finite 
value  of  p  and  disappears  withp-*oo. 

Again  as  in  section  1,  the  following  expressions  are  Introduced: 


X  =  t/t 

S  =  -  1 

Q  =  A  a 

Q  -  A  P  ,  and  furthermore 
o  o' 

_ 

^  “  ef  p  “  fi  P 
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This  gives 


e.  ^ 

Q(1  +  Sx)  -  Q^(l  -  x)  -  A^tix  e'^^j  oe^  dx  =  0 


(60) 


which  for  p -►oo  (  t)  ■♦0)  approaches  equation  (8). 

Now  the  integral  equation  is  transformed  into  a  differei-tial  equa¬ 
tion  by  dividing  It  by  x  e“^*  ,  differentiating  it  with  'respect  to  x, 
and  multiplying  it  by  x^  e~^^,  which  results  in 


^  (X  +  Sx^)  -  (Q  -  Q^)(l  -TlX  +ir)X^)  =  0. 

do 

Since  Q(x)  and  ^  must  approach  the  formulas  (9)  and  (10)  for  x  s  0 
(because  in  the  loginning  the  shunt  resistor  is  r  =  0),  we  have  the  fol¬ 
lowing  initial  conditions: 

Q  =  Q^,  and  ^  +  S)  for  x  =  0  . 

From  the  proper  solution  Q(x)  ,  the  normal 1  zed  current  y(x)  is  de¬ 
rived  as  prescribed  by  equation  (6).  Skipping  some  intermediate  calcu¬ 
lations,  the  final  solutions  to  equation  (61)  are: 


for  8=0 


Q  =  Q. 


1  -  X  e 


tjx 


<1  -  |)] 


y  =  (1  -tjx  tTlx^)  “  2^  , 


for  S  ^  0 


Q  =  Q- 


23*  .  3?±g!g 

1  -  (1  t  S)  X  e®  (1  +  Sx)  ® 


2S^tT 


^  X  - 

y  =  <1  +  s)(l  -T)x  +fjx^)  e  (1  +  Sx)  ® 


(62a> 

(63a) 

(62b) 

(63b) 


lor  p-^oQ  (t)e0),  these  equations  approach  the  former  ones,  namely 

(62a),  (62b)  (9) 

(63a),  (63b><»(10>  . 

The  total  charge  output  within  the  normal  pulse  length  x  =  1  is 
*  Q(l)}  or  relative  to  the  maximum  output  : 
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Examples  are  calculated  for  S  =  0  (c^  =  S  =  1  (e^  =  2€^)f  and 

S  =  2(t^  =  3e^).  Since  the  calculated  values  depend  on  the  assumption 
about  ,  the  dielectric  constant  in  the  uncompressed  material,  its  two 

limiting  values  (see  fig,  2)  €,  =  a  =  1000  e  (for  small  field)  and  = 

€  -  50UU  e  (for  maximum  fie  la)  are  both  use8,  the  correct  value  lies  in 

m  o  ’ 

between  a  and  c  .  The  results  are  shown  in  the  following  table.  In 

most  cases ,  the"*  values  based  on  =  a  or  =  €  are  approximately  the 

i  1  m 

same  for  a  given  p  ,  also  the  Influence  of  S  is  sgall.  In  general,  the 
charge  output  is  higher  than  95  percent  if  p  >  10  ‘  ohm  cm,  and  higher 

than  90  percent  if  p  >  3  •  10**  ohm  cm*. 


Figures  12  and  13  show  some  curves  of  the  normalized  current  output 
y(x)  calculated  for  certain  values  of  p  .  Figure  12  shows  the  case 
=  tj,  according  to  equation  (63a),  figure  13  the  case  =  26^,  ac¬ 
cording  to  equation  (63b)- 


7 .  CURRENT  OUTPUT  UNDER  AN  OBLIQUE  IMPACT 

In  an  actual  experiment  the  projectile  may  hit  the  element  with  a 
larger  or  smaller  angle  of  Impact.  In  this  section  the  deformation  of 


f  Recent  experiments  with  barium  tltanate  ceramic  (ref  2)  gave  values 

Q  -  Q(l) 

-  -  =  91  ±  5  ^  for  110  kilobar,  thus  indicating  a  resistivity  of 

% 

5 

1 0  ohm  cm  . 
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Table  1.  Charge  Output  in  Percent  of  the  Initial  Charge,  if  There  is 
a  Resistivity  p  Behind  the  Shock  Front  (for  Short  Circuit 
R  =  L  =  0). 


Q  - 

Q(l) 

percent  for 

p 

- ^ 

In 

) 

o 

ohm  cm 

"i  = 

"i  =« 

=  c 

'f 

=  e 

m 

"i  = 

=  a 

=  ®^f 

= 

~  a 

=  ®Cf 

lo'^ 

100.0 

100.0 

100.0 

100.0  i 

1 

99.9 

100.0 

3-10® 

99.9 

100.0 

99.9 

100.0 

99.7 

99.9 

10® 

1 

99.6 

99.9 

99.7 

99.9 

99.2 

99.8 

3 '10® 

98.7 

99.7 

99.0 

99.6 

97.5 

99.5 

10® 

96.1 

99.2 

97.0 

98.8 

92.6 

98.5 

4 

3*10 

87.7 

97.4 

90.3 

96.0 

77.4 

95.0 

10^ 

67.4 

92.4 

73.7 

88.5 

46.4 

85.6 

3-10® 

26.7 

76.8 

36.2 

66.2 

7.7 

59.9 

10® 

L.9 

45.4 

4.2 

29.5 

0.1 

21.5 

3*10® 

H 

■■  ' 

1.7 

-  ■  ■' 

0.6 
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a  roHisttvi 


the  pulse  due  to  this  angle  will  be  calculated.  The  short-circuit 
assui.iptlon  (R  =  L  =  0)  and  the  linear  fit  to  the  hysteresis  curve 
(equation  (3))  are  used. 


4 


2r 


0 


+ 


Figure  14.  Oblique  shock  front  through  the  element  (angle  of 
shock  front  exaggerated). 


The  shock  front  originates  at  the  edge  of  the  element  that  Is  hit 
first,  and  the  time  is  counted  from  that  moment  on.  The  back  sur¬ 
face  of  the  element  along  that  edge  is  reached  at  the  time  t  *  -^/u  , 
again  called  the  "normal  pulse  length."  Other  regions  of  the%lement 
experience  the  shock  front  with  a  time  delay,  the  opposite  edge  has 
the  maximum  delay  t^;  along  that  edge  the  shock  front  naches  the  back 
surface  at  the  time  tg  +  t^.  The  time  delay  t_  is  connected  with  the 
angle  of  impact  (p  by  the  formula 


ta  =  —  tan  (65) 

with  2r  =  diameter  of  the  element  and  v  =  Impact  velocity  of  the  pro¬ 
jectile.  For  example  with  2r  =  5/8  in.  and  v  =  1400  ft/sec,  an  angle 
of  impact  as  small  as  1/4  deg  results  in  a  time  delay  t^  ■  0.16  psec, 
which  must  be  compared  with  t^  =0.70  pspc.  The  angle  of  the  shock 
front  in  the  material  with  respect  to  the  planes  of  the  electrodes  is 
about  2.5  deg  in  this  example,  since  u^  is  about  15,000  ft/sec. 

In  order  to  start  the  calculation,  the  element  is  thought  to  be 
split  into  small  elements,  each  with  the  area  dA(0  and  the  thickness  6  . 
Each  element  undergoes  the  same  transition  as  the  shock  front  travels 
through  it,  delayed  in  time  however  as  we  go  from  one  edge  to  the  other. 
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Ai‘  -iM'  k  !  riMit  the  a  I i|^  i  C  13  l4  3 £  'ii  r*  Li  t  X  i  i  t  If  >  t  y  s illtl  t  1 

(c  ,  J  5  dit;,  as  stated  abnve)  so  that,  no  rorrcctlon  Is  needed  lor  the 

I  ;k  '  'hat  the  direct.ion  of  travel  and  the  normal  l.o  the-  shock  front  are 
III  '  in  fiarallcl  Besides  being  the  simples',  possible  formula,  the  short - 
<  1  re  LI  1  t  a.ssumption  is  also  very  convenient  hero  since  it  requires  that; 
the  voliagi  drip  i(  ross  the  electrodes  be  zero.  All  smalt  elements  can 
h<  I  e  feiro  be  'reated  as  independent  elements.  The  following  calculations 
are  confined  to  the  case  t  <  t  .  'I’he  current  o'jt.put  pulse  consists  of 
Mini-  i>arls  Tin-  rise  tinie'^(A)  ?rom  t  =  0  to  t  =  t  is  the  time  during 
whii.ti  only  a  part  (increasing  with  time)  of  the  front  electrod*  is  in  con  ■ 
'act  with  the  projectile  and  therefore-  only  a  part  of  the  oliotrent  is  in- 
vi.lviu,  Liurirg  r.he  mam  pulse  (B)  from  t.^  to  ^  both  electrodes  are  fully 
involved.  The  decay  time  (C/  of  the  pul.se  from  to  t^^  f  Is  explained 
liV  'he  situaMnn  that  the  back  surfa-jc  has  boon  reached  by  the  shock  front 
and  only  a  decreasing  area  of  the  back  elect-ode  is  producing  free  charge. 
It  IS  supposed  that  reflections  from  the  back  end  of  the  element  or  from 
I  he  supporting  rod  do  not  disturb  the  process  The  calculation  starts  with 
a  combination  of  equations  (1),  (2),  and  (3i 


I 


«)  +  f-  ?  -  0. 

f 


(66) 


Tlie  distance  ^  of  the  shock  front  from  the  front  electrode  (fig.  14)  now 
depends  on  time  t  and  position 


2r 


- 

3t 


(67) 


All  variables  are  now  functions  both  of  t  and  .  The  solution  of  equa¬ 
tion  (66)  must  be  integrated  across  the  appropriate  surface  elements 
dAO^) ,  which  gives  the  charge  Q(t). 


If  the  surface  of  each  electrode  is  a  circular  area,  we  have 
_ ^2r 

dA()J)  --  2  J-lTX,  -  1%  ,  I  dA  =  Trr^  -.=  A, 

s  id  the  current  is 


I(t)  =  [ dA. 
J  St 


/  f*<y  \ 


(69) 


The  usual  transformations  are  Introduced; 

u  =  ^’/t 
s  o 

S  .  -  1 

X  X  t/t^ 

t  t 

y(x)  =  -  '  I  (x)  =  •-  1  (x) 

o 
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tk  re  I  h('  s'liHioii  cai'  be  after  ‘  rans  iirma  t  .  ops 

H  •  -  X ,  fJi:;  =  -  dx 

>  •> 

2r^  -  a  ‘  bx  ••  x“ 

a  ■  B(,!r  in 

Ii  -  2(fi-r) 

T'tie  solutions  of  the  Intergral  are  different  for  a  >  0 ,  a  0 ,  and 
a  V  0.  Our  calculations  are  confined  to  a  <  0,  which  is  equivalent 

I  () 


1 

X  >  z  •  -  . 

Since  actually  S  falls  between  1  and  2,  this  condition  Is  always  ful¬ 
filled  for  z  <0.5  For  z  >0.5,  the  very  first  part  of  the  pulse 
may  require  a  different  solution.  Here  we  proceed  with  a  <  0,. 

The  general  solution  of  the  integral  in  equation  (72)  is 


T  sin 


.  i 


with  the  abbreviation 


^  Bin'll- 

(73) 


^  _  _ l4^Sx-Sz/2 _ 

■«/B(B-2r)  s/d+SxXl+Sx-Sz) 

710  complete  the  solutions  (71)  or  (73),  the  limits  cf 


(74) 


the  integrals 


/ 


^^2 


are  to  be  put  in.  These  limits  depend  on  the  region  of  the  pulse,  s.z 
explained  above  in  this  section. 


Region  A  (rise  time,  0  <  t  <  t^,  0  <  x  <  z):  for  a  fixed  time  t 

(or  x)  within  this  region,  the  front  electrode  is  covered 

„  j  2r  2r 

by  the  projectile  between^  1  "  ®  ^2  '  t  ^  *  z  * 

a 

Region  B  (main  pulse,  t  <t<t  ,  z<Xvl);  both  electrodes 
are  fully  involved,  so  =  C  arc*  _  2r. 
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Region  C  (decay  time,  t  <t<t  +t,  l<x<l+z):  the 

back  electrode  is  only  partly  involved,  therefore 

^  ,  =  J-  (L  -  t  )  =  —  (X  -  1),  =  2r  . 

1  t„  o  z  z 

a 

Now  the  current  output  for  the  three  regions  can  be  written. 

For  S  =  0: 


4x  -  2z  , — ; - 1  -l  x 

- ^  ^/x(z  -  +  -  sin  (2  - 


1'  ^ 


for  0  <  X  <  z 


(71a) 


y<x)  =  1 


for  z  <  X  <  1 


(71b; 


y(x)  =  -  -  i)(z  -  X  1) 


(71c> 


-  -  sin‘^  (2  —  -l)^iforl<x<ltz 

TT  Z  2  — 

Figure  15  demonstrates  y(x,  S«0)  for  z  s  0  and  z  s  0.4. 

For  S  ^  0  (using  (73)  and  (74)): 

“  x)  (1  -  2  |)  ♦  ^  (T  -  1) 

IT  S  z 


-  T  sin 


1  -  2  I  (1  -  Sz  +  Sx)|  for  0  <  X  <  z  (72a 


y(x)  »  (T  -  1) 

S  z 


for  z  <  X  <  I  (V^v*) 


y(x)  .  -  ^^/(x  “  l)(z  -  X  +T)  -  sin‘^  (1-2 

IT  z^  L 

♦  T  sin"^  (^  “  2  5^  |l  -  +  I  ‘  i>] 

for  1  <  X  <  1  +  z  . 

(72c: 
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Figure  16,  Influence  of  ttie  angle  of  impact-  (z-'-tan  c|7)  on  The  curve  shape 
(for  €  ”2fR=L=0>. 


Figure  16  demonstrates  y(x^ 


S  =  1)  for  z  =  0,  0.2_,  0.4,  0.6  0  8, 

1.0, 

8.  THE  SHORT-CIRCUIT  ASSUMPTION 

The  mathematical  treatment  of  many  features  (such  as  the  curved 
hysteresis,  maximum  energy  oulput,  conductivity  behind  the  shock  front, 
oblique  impact)  is  greatly  facilitated  or  even  only  possible,  if  the 
short-circuit  assumption  R  =  L  =  0  is  made.  Our  experiments  use  an 
electric  load  R  =  10  ohms  and  L  <  O.S  •  10““^  henry,  and  the  question 
has  to  be  investigated  whether  the  short-circuit  formulas  can  bo  used 
with  these  values. 

Figures  4  through  7  show  how  close  curves  with  these  values  of  R 

and  L  come  to  the  limiting  case  R  »  L  =  0.  Oniy  the  initial  part  of 

the  pulse,  for  0  <  t  <  0.1  •  t  ,  is  affected  at  all. 

—  —  o 

Another  way  of  justifying  the  short  circuit  assumption  with  R  =  10 
ohms  is  to  look  at  the  element  as  a  capacitor  (capacitance  C  =  7“10~^^ 
farad).  With  a  short  circuit  the  charge  flows  off  iaaediacely  after  it 
is  set  free  (RC-time  *  0).  In  our  case  the  ROtime  is  7- 1.0“®  sec,  which 
is  only  1  percent  of  the  normal  pulse  duration  t^,  *  7*10'‘'^  seconds. 

A  third  way  of  showing  that  R  s  10  ohms  is  small  enough  to  Justify 

the  short-circuit  assumption  is  to  use  equation  (1>.  As  the  curves  for 

small  R  show  (fig.  5),  for  a  time  x  *  0.3,  (traveled  distance  5  =  0.1  cm) 
about  50  percent  of  the  total  charge  has  flowed  off  the  electrodes  (o(O.d) 
s  Pq/2).  At  this  moment  the  voltage  between  the  shock  front  and  the  back 
electrode  is  approximately,  according  to  equation  (1): 


1  ^i 


0.5- (8) -IQ-* 
0.886*10“^® 


(0.22) 


-  9900  V 


At  the  same  time  the  voltage  between  the  shock  front  and  the  front  elec¬ 
trode  is  approximately,  according  to  equationd) : 


0.5P 

o 


O.S‘(8)-10~® 

0.443-10“^® 


(0.1)  s  9100  V 


The  current  at  this  time  is  only  about  25  amp  (measured  value),  and  the 
voltage  across  R  therezui-e  250  v.  This  demonstrates  that  neglecting  the 
third  term  in  equation  (1)  does  not  introduce  a  big  error. 
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METHODS  FOR  EVALUATING  Gf  FROM  EXPERIMEIfTAL  DATA 


The  equations  in  the  previous  sections  contain  all  or  some  of  these 
values;  A,  t,  P  ,  Q,  t  ,  R,  L,  p ,  a,  3, <P ,  ,  and  e^.  If  Q  approaches 

A‘P  ,  then  the  resistivity  P  is  sufficiently  high,  as  shown  in  section  6, 
and°we  need  not  correct  for  any  current  leakage  through  the  compressed 
region.  Then  everything  can  be  measured  before  or  during  the  experiment, 
except  Gf.  This  quantity  therefore  is  determined  by  the  outcome  of  the 
experiment,  i,e.,  the  current  output  y(x)  or  I(t).  is  the  dielectric 
constant  of  the  material  under  high  compression  and  elevated  temperature, 
the  temperature  being  determined  along  with  the  compression  by  the  methods 
discussed  in  reference  1.  Cf  is  an  important  property  of  the  material 
under  those  conditions,  and  it  is  one  of  the  main  objects  of  our  present 
experiments  to  get  values  for  ej .  Most  of  the  formulas  of  the  preceding 
sections  are  too  complicated  to  be  solved  for  Gj,  and  we  have  to  look  for 
the  most  simple  expressions  and  elect  the  appropriate  experimental  arrange 
ments.  One  of  these  is  the  open  circuit  situation 


y(x) 


“  GjAR 


from  which  follows 


U(t)  =  R-I(t) 


"f^o 


and 


dt 


"f^o 


This,  by  the  way,  is  the  initial  rise  of  U(t)  for  any  R  (equation  19) 
which  for  an  open  circuit  goes  on  as  a  linear  rise. 


Equation  (75)  is  very  simple,  and  one  need  only  measure  the 
rise  of  the  curve:  voltage  versus  time.  Hie  difficulties  however  are: 


1)  Very  high  voltages  are  reached  within  a  short  time  (with  our 
experimental  values:  60  kv  at  t  =  t^  =  0.7  ^sec),  which  may  genera've 
breakdowns  through  the  element,  around  the  edges,  and  across  the  probes. 

2)  t  is  not  easily  found,  since  the  actual  scope  traces  of  open- 
circuit  experiments  do  not  show  a  clear  picture  of  the  pulse  and  especially 
its  end  part. 

3)  Disturbances  of  the  rise  time,  due  to  the  angle  of  impact,  to  the 
amplifiers,  or  to  any  other  reason,  are  superiMiH^sed  and  therefore  falsity 
the  resulting  apparent  value  of  €^. 

It  might,  however,  be  possible  to  overcome  these  difficulties  and  to 
get  proper  results  by  this  method. 
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The  ot*  r  simple  arrangement  is  the  short-circuit  situation,  which 
is  treated  in  many  details  In  this  report,  and  which  is  nearly  fulfilled 
in  our  experl  lents  (aee  sec-.iou  8).  The  rest  of  this  section  (9)  deals 
with  method*'  ;:o  eval*- ited  €,  from  a  given  curve  I(t>  or  y(x)  under  short 
circuit  conditions. 

Usually  a  eur’-e  U(t)  is  obtained,  with  a  known  value  of  R  (10  ohms). 
Therefore  I(^)  is  given.  This  curve  generally  shows  a  wain  pulse,  past 
which  the  curve  goes  through  the  sero  line.  Hie  current  then  bt cones 
negative  for  a  while,  followed  by  a  sealler  pulse  and  ripples.  The  pulse, 
as  far  we  need  it,  is  at  its  end  at  the  crossover  with  the  zero  line. 
Everything  that  follows  is  influenced  by  reflections  and  is  of  no  interest. 
The  charge  output  is  Q  dt  across  that  ealn  pulse,  the  experlaental 
values  of  Q  are  close  to  (^,  and  are  put  equal  to  for  the  following  cal¬ 
culations.  The  pulse  never  appears  to  be  ideal,  the  rise  and  decay  are  not 
vertical  but  inclined  lines,  even  where  a  snail  angle  of  inpact  has  been 
measured.  Therefore  the  nonoal  pulse  length  t  has  to  be  taken  between  (a) 
the  half  value  of  the  first  straight  line  of  tSe  rise  and  (b)  the  half 
value  of  the  decay  of  the  pulse.  This  is  in  agreenent  with  figures  8 
through  10,  which  show  that  the  rise  of  the  curves  consists  of  a  first 
straight  line,  followed  by  a  bend  and  a  nore  Inclined  part.  In  figure  17, 
the  Important  parts  of  the  curve  are  identified. 


PEAK 


Figuie  17.  General  pulse  shape  for  short  circuit. 


With  t^  and  Q^,  the  transfomation  to  the  nomallsed  form  x  a  t/t  , 
y  =  ^  made 
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Due  to  the  deviation  from  the  ideal  curve  shape,  as  stated  before, 
the  leading  part  of  the  pulse  is  not  well  suited  to  derive  Cf  from  it. 

The  tail,  however,  follows  closely  the  theoretical  formulas  except  for 
the  decay  part  of  the  pulse  which  also  cannot  be  used.  Tlie  quickest, 
although  not  very  precise,  method  to  evaluate  cy  appears  by  use  of  equation 
(10),  which  solved  for  e^Ay  =  S  +  1  gives  us: 


•  2xy(l-x)  -  -^/l  -  4xy(l-x)  J  /2x^y 


El/€, 


CiAj 


€,/Cj 


=  [•- 


2xy(l-x)j/2x  y 
2xy(l>x)  +*Vl  -  4xy(l>x)j  /2x^y 


for  J  <  X  V  X 

—  Hi 


for  X  =  X 


(76) 


1/y 


for  X  <  X  <  1 


for  X  ■  1 


with  X  m  -i-  ■  — - 

m  2-f8 


is  a  mean  value  of  the  dielectric  constant  in  the  uncompressed  material, 
as  explained  in  section  1. 


A  better  method  is  to  use  equation  (34)  instead  of  (10),  thus  employ¬ 
ing  the  curved  hysteresis  i*«etead  of  the  linear  fit.  Hals  brings  us  nearer 
to  the  actual  pulse  shape.  In  doing  so,  the  quadratic  fit  a,  to  the 
hysteresis  curve  (equation  4)  must  be  known,  instead  of  only  the  mean  value 
Cy.  The  general  formula  (34)  cannot  be  solved  analytically  for  Cy,  an 
approKlmate  formula  must  be  invented.  This  is  done  by  creating  an  auxiliary 
function  <^(x),  which  derives  it  Justification  solely  from  the  fact  that  it 
approaches  the  tail  of  y(x)  so  closely,  as  can  be  seen  in  figures  8,  9,  and 


A  curve  (similar  to  equation  10)  y(x)a 


[l4-(k-l)x 


^  for  0  <  X  <  1, 


is  transformed  to  the  range  x  <  x  <  1  by  the  linear  transformation 

' 


giving 
y(x)  «  ^(x) 


SI 


The  value  of  (j>(l)  =  ^  is  forced  to  coincide  with  y(l>  =  t./e  (from 

f  m 

equation  34).  This  determines  k  =  and  thereby 

.2 


(^(x)  = 


V. 

7~ 

€f  P 


I  Cx 


€m-ef 
X  f  - 

P  e# 


2  for  X  <  X  <  1  , 
P  - 


(77) 


V  is  the  position  of  the  peak  of  the  curve  y{x)  according  to  equa- 

ti'on  (34).  €  is  defined  with  equation  (34).  Equation  (77)  can  be 

solved  for  e  /e.,  giving 
m  I 

[  y  ■  ♦ 


By  putting  the  values  of  x  of  e  and  any  values  of  x  and  y  (the 

latter  for  (^)  from  along  tall”of  the  measured  and  normalized  out¬ 

put  curve  y(x)  into  equation  (78),  this  equation  yields  c*.  Figures 
8,  9,  and  10  show  which  part  of  the  tall  is  suitable.  * 

ACKNOWLEDGMENT 

The  author  wishes  to  thank  P.  S.  Brody  for  introduction  into  this 
field  and  for  helpful  discussions,  and  H.  R.  Kollmeyer  who  did  most  of 
the  work  for  the  experimental  setup  and  helped  carry  out  the  experiments. 
R.  F.  Butler  set  up  the  suialog  computer  which  drew  some  curves  lor 
section  3.  Hie  work  was  supported  as  part  of  a  Transducer  Reseasich 
Project. 

REFERENCES 

1)  P.  S.  Brody,  "strong  Shock  Waves  in  'Polled'  Barium  Titanate  Ceramic 
Elements,"  DOFL  Report  1R-869  (October  20,  1960). 

P.  S.  Brody,  R.  B.  Wlttckxndt,  "Dielectric  Constant  of  Barium  Titanate 
at  100  Kilobars,"  DOFL  Report  TR-917  (March  15,  1961). 

3)  R.  Schall,  K.  Vollrath,  "Zum  Verhaltcn  fei'r  •^lektrischerKeramik 
bel  intens^ven  Stosswellen,"  lumri  .it^.uber  <.tr.  IV,  Internatlonalen 
Kongress  fur  Kurzzeltphotographie,  Xbln  1958,  p  .329, 

4)  E.  Kaake,  "(iAwohnliclie  Oifferentlalglrlchungcn,"  Akadamlsche  Ver- 
lagsgesellschaft,  Leipzig  (1943). 


52 


DISTRIBUTION 


lie  par  line  r.‘  of  the  Ar;r.y 

Office  of  Lhe  Chief  of  O’ dnance 

file  Pentagon,  Washington  25,  D.  C. 

Attn  ORDTN  (Nuclear  &  Special  Components  Branch) 
Attn’  ORDTB,  Joseph  Kaufman 

Director,  Army  Research  Office 

Office  of  the  Chief  of  R?search  &  Development 

Department  of  the  Army 

Washington  ,  D.  C. 

Commanding  Officer 
Ficatlnny  Arsenal 
Dover ,  New  Jersey 

Attn;  Felcman  Research  &  Englreerlng  Laboratorle-:. 
Attn;  Library 

Commanding  Officer 

U.S.  Army  Signal  Research  &  Development  Laboratory 
Fort  Monmouth,  New  Jersey 
Attn:  Tech  Library 

Commanding  General 
Aberdeen  Proving  Ground,  Maryland 
Attn:  BRL—R.  J.  Eichelbt*rger 

Commanding  Officer 

U.S.  Army,  Office  of  Ordnance  Research 
Box  CM,  Duke  Station 
Durham,  North  Carolina 

Commanding  General 

Engineer  Research  &  Development  Laboraloiles 
i:  S.  Army 

Fort  Beivcir,  Virginia 

Attn:  Tech  Documents  Center 

Commander 

U.S.  Nava.  Ordnance  Laboratory 
Corona,  Ca-ifornia 

Attn:  Documents  Librarian 

Commander 

U.S,  Naval  Ordnance  Laboratory 
White  Oak,  Silver  Spring  19,  Maryland 
Attn;  M.  Solow 
Attn  S.  J.  Jacobs 
Attn;  B.  Drimmer 


DISTRIBUTION  (Continued) 


Depo.rtmc-nt  of  the  Navy 
Bureau  of  Naval  Weapons 
Washington  25,  D.  C- 

Attn:  DLI-3— I’ech  Library 

Commander 

Naval  Research  Laboratory 
fYaS  hington  25,  P.  C. 

Attn;  G.  Irwin 
Attn:  Tech  Library 

Commande  r 

Armed  Services  Technical  Information  Agency 
Arlington  Hall  Station 
Arlington  12,  Virginia 

Attn:  TlPDR  (10  copies) 

Sandia  Ckirporatlon 
Sandia  Base 

Albuquerque,  New  Mexico 
Attn;  G.  Anderson 
Attn:  F.  Nielson 
Attn:  R.  Grahm 

Attn:  H.  G.  Baerwald 

Los  Alamos  Scientific  Laboratory 
Lub  Alaittos,  New  Mexico 
Attn:  J.  Wackerlee 
Attn:  J.  M.  Walsh 
Attn:  R.  Duff 
Attn*  S.  Minshall 

University  of  California 
Radiation  Laboratory 
Llverqtore,  California 

Attn:  Taylor  Abegge 

Director,  Stanford  Research  Institute 
Menlo  Park,  California 

Attn:  The  Poulter  Laboratories 

Clevite  Research  Corp 
Cleveland,  (Hiio 

Attn:  H.  Jaffe 
Attn:  D.  Berlincourt 

National  Bureau  of  Standards 
Mineral  Products  Division 
Industrial  Bldg. 

Attn:  A.  D.  Franklin 


54 


DISTniCUTIOX  <Cor.tir.ued) 


National  liureau  of  Standards 
"hoi'tnody  tiartii  cs  Di  vl  si  on 
Building 

All.  n:  M  S  Orfen 

National  Bureau  of  Standards 
Mechanics  Division 
Industrial  Building 

Attn:  E,  C.  Lloyd 
Attn;  L.  K.  I rwin 


Internal 

Hinman,  W.  S.,  .Ir./McEvoy^  R.  W. 

Apsteia,  M./Gerwin,  H.  L./Guarino,  P.  A./Kalmus,  H.  P. 
Sures,  A,  II  /Schwenk,  C.  C. 

Hardin,  c.  D, ,  Lab  100 
Horton,  B.  .M  ,  Lab  200 
Rot kin,  I . ,  Lab  300 

Landis,  P.  E ./Tucci nardi ,  T.  E.,  Lab  400 
Hatcher,  R.  D. ,  Lab  500 
Flyer,  I.  N. ,  Lab  fiOO 

Campagna,  J.  H./Apolenia,  C.  J.,  Dlv  700 
DeMasl,  R.,  Dlv  800 

Franklin,  P.  J. /Horsey,  E.  F.,  Lab  900 

Harris,  F.  T. ,  310 

Bowles,  R.  E.,  310 

Brody,  P.  S.,  310 

Kinzelman,  G.  W  ,  320 

Piper,  W.,  320 

Ravitsky,  C. ,  310 

Warren,  R.,  310 

Curchack,  H.,  310 

Horn,  L.,  310 

Wlttekindt.  R.,  310  (20  copies) 

Seaton,  J.  W. ,  260 

Technical  Reports  Unit,  800  (3  copies) 

DOFL  Library  (5  copies) 

Technical  Information  Office,  010  (10  copies) 


(Two  pages  of  abstract  cards  follow.) 


I 


K5  ra 

SI  2  3  2 


9  S!  e 

9  a  si 


x  !•  a  ' 

8  a  &i5 

£  0  O  ■ 

m  o  XT  < 

M 

«  M  «l 

•*  V 

>,  <%  <c 

'  8.  §  ♦' : 

T>  fV  »<  «» 

4)  « 

«  4> 

fc*  ♦«  «  £ 

>  &  M 

tt  **  U  0 
%•  a  0  91 

t  0  P  0 


m  B.  Ji  ^ 
4*  C  H 

^  2  .  P 

0  9  0  P 

e  0  o  0 

4>  M  ^  0 

^  H  c  n 

•  .8  E  u 

a  c  »«  0  c 
—  0  ?  •*  '■ 
h  -4  0 
**  4*  ^ 

U  <4  *>  <•  « 

*»  H  -2  I 

<H  iS  •  • 

tt  fc«  V  C  ( 

5  «  0  1 

b  ^  ^  M 

Coo  I 

•  &  ‘0  4 

•to  •  C  ' 

•»  >  «  4 

r  4J  C  _  ' 
N  0  •  < 

«4  «  M  ' 

b  bi  «  «4  ' 

s  zi-i 

.3*8 

a  •<  &  4 

i 


as  ra 

ss  2  3  2 

M  ^  »to 

B  O  #4  M 

-•S  Sd 


J«  O  3  < 

4#  #5  »• 

o  K*  n  y 
f  o  o  • 

M  O  «  • 

M 

-  rs .: 

X  ^  £ 

‘’88-; 
■o  n.  to  o 

4)  «  «l  b 
M  O 

b  b  «t  x: 

O  3  3  i» 

>  a  «i 

«  b  b  o 


.to  »to  b  3  O  •«  M 

#4  M  ««  O  >  4  O 

M  V  *• 


•  •  4 

»  3  b 
I  O’  b  3 
^  to  3  U  C  * 

:  ^  o.  « 

A  M  3 

i  o  g  o^. 

*6  b  b 

4  X  e  • 

>  P  *2 

:  «to  b  b  3 

I  0  *  b  u 

r  3  ^ 

K  *3  •  U  «• 


»  •  13  ♦*  b 

1*3  S  b  ^ 


I  ^ 

9  a 

lwP«V>» 


4  b  b  4 
b  e  3  0  6 
•  aA- 

4i  b  b  b 

e  b  a  0 
b  3  o  b  £ 

fi  b  Off 

b  b  b«  G 
b  •  O  « 

•  •  f- 

to  b  b  G 

u  e  b  o  c 
b  o  3  ^  ^ 
b  b  U 
b  b  b  *• 
U  4  «  4  « 

•  N  b  b  I 
b  b  C  •  1 
4  b  b  C  ( 
0  4  VI 
b  b  b  tj 

b  Q  O  I 

0  O.  <t  I 

b  OR' 

•  >  4  « 

3  A  b 

•  b  3  3  ( 
M  0  V  ’ 
b  4  b 
b  X  O  b  * 

4  p  A  fb 
b  b  b  4 

is*.  E 
-3  * 

R  M  &  4 
0  4  4 

.R  x:  js.  u 

;•  M  VI  -to 


;  f.  w  to  to  to 

I  b  4  O  4  0  %• 

)  b  S  v«  o  X 

o  b  b 

I  4  O  R  e  t  T* 

•  b  X  O  4  A  b 

b  K  b  • 

I  3  e  .«  •  3*  b 

f  O  O  4  0  4 

4  M  to  9  *44 

•  4  b  b  4  9 

4  9  b  (b  4  3  3 

»  b  3  y  0  o 
n.  R  4  b  •  4 

s  i!  o  *'2  2 

:  3  “SS  *  £ 

O  O  4  4  a 

0  A  P  b  U  ^  ^ 

b  4  AO 

to  r.  *1  b 

0  4  *0  3  0  0 

i«  b  3  A  O  Si  A 

0  O  b  b  b 
to  X  4  4  b 


I  b  4  •  O  I 

I  b  b  3  A 
»  O  3  O  b 
)  3  4 

I  3  4  4  b 

I  4  b  3  O 

O  O  V« 

I  O  b  4 

I  O  b  4 

1  4  4  3  O 

.  a  > 


0  4  114 

b  V  9  U  M  9 

X  9  b  9  •  ^ 


X  b  •  2  •  ^  ^ 

51  ia.  ‘IS  1-5^2- 
|P I 8S“S3 i-3’ i 
•  .®l*'o.8s8toia 

*i-a-8“S>.ss-fl 

3&2p4S40|MObJ 
oabbgM  ._4S!^i2 


*  a  i 

8  d  a 


?  &  *  i8*'S. 

23 

toO  bbgAOlB 


1 25  *1  It  ilUll^Sis 


? 

4 

•  t 


i8  5| 

25  M 

S  w  X 

H  *4  *4 

6?  -2 
X  M  e 
U  &  M  91 


bbXM  bto  f|4b 

S&*«S..3£dl831 

<S  •4bj»3  i^b  b 

vJKC^uS'^bSjlg  *4 

2-z:28  a^se  ^2- 


lx  fl 

o  c 

B  -<4  C 


p  c  I  I  •  2  ft  • 

E  'OOUMC  xt  V 

M  * 

« six • «: 

r§H:-3n-lJ 

;  ;  >  O  ®  2  “  3  *  ?  5  8 

*&8“::S2~sCa35 
l‘}'S233“  isi-a3 

S''*'5cC«TiS6a3x 

5l2*2225|i;|l 

-**83^|53|-&S 

«I4«M  4  '*4 

-  -  ia2 


gi 


u  rt  «« 


:g 


2  IS  i  « 

"  ^8  it  I'hrisii'li- 

!S2=|Hs:J5i«lJEi 


o'  ' 
9  o 


C  I  I  <A 


•f  V  f* 

Q  tt  9/  *  ** 


:  «4  u  u  3  «« 


9  U. 


s  Ci  t.  **  ft  xi  O  C 

UJ  p  9  9  if  C  ^ 

I  M  >  O.  « 


*«  k  K  O 


U  ^ 
x4  l«  O 


«*  •>  n  u  o 

MO  ftt  W  Vl 


-2 


«•  M 


iCtO 

g  C  O  *•  **  £  U 

I  C  V  C  *  V  9  9  # 

»e-«go  •* 

>  >.Mj^la 


a 


<  « 

f-  U 
vt  o 


It  M  o  «* 


>  5  ?  t 

«  .«  £  It  9 

^  H  •  ••  0 


£ 


*4  ^ 


S  rfi 


B*  ^ : 

«f  «•  I 


9  9.  « 

^  M  M 


3 1  i 


::B  "lii 
B  s  E^- 
-.83 

BO  •  ^ 


88S,|5 

2  •  V  M  9  «  ; 

JM  K ««  2  O  e  ' 

o  *a*3  *? 
**5-8  *.85 

is 

82 • • 8  52 
!S52a*®8; 

I  a  <»4  0  9  ■  t  _  ' 

9,  •  5  «  2  M  C  ' 

9  P  S  9  M  9 

t  ft  r*  3  ®  * 

«  0  It  i  *t 

llEaai .! 

I  C  «f  «*  ««  B 

:&8888p 

I  *  5  2  :  2  Li 


!l 


I  p  «i  *a!?|i«»;sl 

<pio  .<x*«Cx«< 

-5>^  a*#*:"? 


9  Z 


9  0  •*  9  ^ 

M  ca  «t  *»  9  9  • 

i:iS5i:ih 

2C*Upa  xaflx 

-8888256  *1 

•  Nxgkxx'fi 

Sas*  *55582 
5-a"82|«:* 

P'lli'i’! 

i«f9«t9«f**0««««li 


